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Q ■ ABSTRACT 

C/5 ' In this work, starting from 21 months of data from the Fermi-Large Area Telescope 

03 ; (LAT), we derive maps of the residual isotropic y-ray emission, a relevant fraction of which 

is expected to be contributed by the extragalactic diffuse y-ray background (EGB). We search 
' for the auto-correlation signals in the above y-ray maps and for the cross-correlation signal 

^ I with the angular distribution of different classes of objects that trace the large scale structure 

of the Universe. We compute the angular two-point auto-correlation function of the residual 
\G ' Fermi-LAT maps at energies £ > 1 GeV, £ > 3 GeV and £ > 30 GeV well above the Galac- 

OO , tic plane and find no significant correlation signal. This is, indeed, what is expected if the 

'nI" ■ EGB were contributed by BL Lacertae (BLLacs), Flat Spectrum Radio Quasars (FSRQs) or 

I star-forming galaxies, since, in this case, the predicted signal is very weak. Then, we search 

. for the Integrated Sachs-Wolfe (ISW) signature by cross-correlating the Fermi-LAT maps 

with the WMAP7-Cosmic Microwave Background map. We find a cross-correlation consis- 
tent with zero, even though the expected signal is larger than that of the EGB auto-correlation. 
Finally, in an attempt to constrain the nature of the y-ray background we cross-correlate the 
Fermi-LAT maps with the angular distributions of objects that may contribute to the EGB: 
^ ! QSOs in the SDSS -DR6 (Sloan Digital Sky Survey- Data Release 6) catalog, NRAO VLA 

H ' Sky Survey (NVSS) galaxies. Two Micron All Sky Survey (2MASS) galaxies and Luminous 

. . .' Red Galaxies (LRGs) in the SDSS catalog. The cross-correlation is always consistent with 

zero, in agreement with theoretical expectations, but we find (with low statistical significance) 
some interesting features that may indicate that some specific classes of objects contribute to 
the EGB. A analysis confirms that the correlation properties of the 21 -month data do not 
provide strong constraints of the EGB origin. However, the results suggest that the situation 
will significantly improve with the 5- and 10-year Fermi-LAT data. The future EGB analysis 
will then allow placing significant constraints on the nature of the EGB and might provide in 
addition a detection of the ISW signal. 

Key words: cosmology: theory - cosmology: observations - cosmology: large scale structure 
of the universe - gamma rays: diffuse backgrounds 



1 INTRODUCTION 

One of the most interesting open problems in astrophysics is 
the nature of the diffuse y-ray background at GeV energies. The 
presence of a diffuse signal was first detected by the OSO-3 
satellite l lKraushaar et al.l 1972h with indi cations of both Ga lactic 
and isotropic diffuse emissions. SAS-2 ( lFichteletal.lll973h and 



later EGRET jSreekumar et al.|[l998h revealed more clearly the 
isotropic component, which is commonly known as extragalactic 
y-ray background (EGB). Although several l ocal processe s have 



been proposed to explain thi s background (iKeshet et al, 



2004 



iMoskalenko et al.l |2006|. |2007| : iMoskalenko & Ported \2Q0% . the 



EGB is generally believed to be the superposition of contributions 
from unresolved extragalactic sources and, perhaps, diffuse GeV 
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emis sion processes lik e the annihilation of dark matter (DM) parti- 
cles dUllio et alj|2002l) . We will, however, in the following use the 
abbreviation "EGB", even if the extragalactic origin of this compo- 
nent is still not fully clear. 

Blazars are the most numerous population detected by 
EGRET and the Fermi Large Area Telescope (Fermi-LAT) 
( iHartman et ar]|l999l : IXbdo et al.ll2009l) . They have long been re- 
garded as the most likely candidates to provide the bulk of the 
EGB emission. However, the recent analysis of the First Fermi- 
LAT AGN catalog obtained after '-^ 1 year of data taking revealed 
that blazars can contribute 23 ± 5%(stat. ) ±12%(syst.) of th e EGB 
in the range between 0.1 and 100 GeV jAbdo et al]|2010dh . disfa- 
voring the interpretation of the EGB as mainly consisting of un- 
resolved blazars. DM candidates in supersymmetric theories (as 
well as other WIMP candid ates) can annihilate into GeV photons 
and contri bute to t he EGB (UUio et al .1120021 : lAndo et aljT2007bl : 
IPieri et al. 2009; Za vala et al.l2010l . l201 l b. However, the amp htude 
and shape of the observed EGB spectrum jAbdo et al]|2010al l3) to- 
gether with the available multi-wavelengt h and multi-messenger 
astrophysical constraints jPato et al.ll2009h seems to indicate that 
y-ray photons from annihilation of WIMP-like Cold Dark Matter 
particles provide, at most, a minor contribution to the EGB. Finally, 
the EGB could be mainly consisting of a population of numerous 
but faint y-ray sources: normal, star-forming galaxies are typical 
candidates that could make a substantial contribution to t he EGB 
below 10 GeV l lPavlidou & Fieldsll2002l : [Fields et al.ll201(]|) . 

To tackle the problem of the EGB nat ure one can use differ - 
ent and rather comple mentary methods (e.g. lHenslevetal1 ( l2009l) : 
iDodelson et"ai] ^20091) ) that may serve as independent constraints: 

• Resolve the largest possible fraction of the EGB into individ- 
ual sources. If the EGB were mainly contributed by a population 
of rare, bright objects, then Fermi-LAT will eventually be able to 
resolve a significant fraction of this radiation and to disentangle 
possible multiple components. If the EGB were mainly contributed 
by a population of common, faint object s, then only a few of them 
will be resolved by Fermi-LAT (see e.g. IPavlidou & Fieldsl ( 1200 ih 
for the case of star-forming galaxies). Finally, if the EGB were sig- 
nificantly contributed by DM annihilations, then, the Fermi-LAT 
might be able to resolve only some individual sources associated to 
foreground dark Ga lactic s ubhalos, depending on the precise nature 
of the DM particle jPieri e t al. 2008, 2009). 

• Compare the observed EGB energy spectrum with model pre- 
dictions based on the luminos ity function of som e classes of objects 
and their energy spectra (e.g. [Fields et al.lfcOlOh ). 

• Analyze the one- or two-point statistics of the observed EGB 
photon counts to distinguish the contribution of clumpy compo- 
nents, typic ally associated to individual sources, to that of a diffuse 
component ()Lee et alj|2009l : [Fields et alj[2010l : [Baxter et al.[[2010l : 
[Slatver & FinkbeineJ[2oioh . 

• Analyze the angular correlation properties of the EGB and 
compare it with those of a pop ulation of y-ray em itting objects 
like blazars ([Ando et al. [2007bl la[). galaxy clus ters jMiniati et al 
20071 : [Ando et alj [2007bl) ty pe la su pemovae teh ang & Beaconj 
2004) star-forming galaxies (I Ando & Pavlidou 200*^ and DM ha- 



los jAndo & Komatsu"2006': Ando et al. 2007b: Cuoco et al. 2007 



Hooper & Serpico . 2007: Cuoco et al. 2008: Siegal-Gaskins, 20081: 
Taoso et al.[ [2009i: iFomasa et al.[ [2009[: [Andd l2009l: [Cuoco et^ 



20101) . 



to the presence of a well defined population of objects. Due to its 
isotropy, the ACE signal of the EGB is expected to be weak. The 
signal can be amplified by cross-correlating the EGB with the an- 
gular distribution of the supposed EGB sources, since the number 
of contributing sources is typically large. For this reason we cross- 
correlate the EGB with different catalogs of extragalactic objects 

that can contribute to the EGB. 

Our approach is very similar to that of [Ando & Pavlidoul 

( l2009l) but our analysis is more comprehensive in two aspects. 
First of all, since our aim is to identify the contributors to the 
EGB among a variety of possible candidates we cross-correlate 
the EG B maps with several d ifferent objects' catalogs (SDSS-DR6 
quasars (iRichards et al.l2009h . SDSS-DR6 Luminous Red Galaxies 
jAbdalla et al.[[2008[) , NVSS radio galaxies ( [Blake & Wal"i|[2002[)) 
in addition to the 2MASS star-forming galaxies djarrett et alj2000l) 
considered by [Ando & Pavlidoul ( l2009h . Second of all, not only we 
do provide theoretical predictions of the expected ACE and two- 
point cross-coiTelation function (CCF), but we actually estimate 
them from the 21 -month Fermi-LAT data. Finally, we perform our 
analysis in configuration space (although we also show the angu- 
lar power spectra for compl eteness ) and consider energies £ > 1 
GeV, slightly higher than in Ando & P avlidou (20oJ), to reduce the 
contamination of the signal due to the Galactic foreground. 

In our work we also pursue another important goal: we attempt 
to det ect the Integrated Sacks-Wolfe (ISW) signal l lSachs & Wolfd 
[1967'). This effect is related to the variation over time of the grav- 
itational potential <I> which arises at late cosmological times when 
when dark energy or curvature become important, and induces ad- 
ditional anisotropies in the Cosmic Microwave Background (CMB) 
at large angular scales. The effect is instead absent during the mat- 
ter dominated era when the gravitational potential is constant, and, 
for this reason, is a potentially powerful probe for dark energy. 
Successful searches for the ISW effect have been performed in 
the past by cross-correlating the CMB maps with the large scale 
stmcture s probes that we have already mentioned (NVSS radi o 
galaxies dNolta et alj2004[:[Raccaiielli et alJ[2008[:[Xia et alj2010l) . 



2MASS galaxies (jAfshordi et al 



2004 



Rassat et alj[2007[). SDSS 



quasars ([Giannantonio et al.l2006[ : [Xia et al.l2009[) ). SDSS galaxies 



In this work we follow the last approach and estimate the two- 
point angular auto-correlation function (ACE) of the Fermi-LAT 
EGB to identify those characteristic features that can be related 



dCabre et al 2006'. '2007^ and with combinations of different trac- 
ers ([Giannantonio et al. 2008: Ho et al. 2008: Hirata et al. 200^). 
Besides the above surveys, also the EGB sources are a probe of the 
large scale structures, and cover a large fraction of the sky. Thus, 
we do expect some ISW signal which could be detected through the 
cross-coiTelation between Fermi and CMB data. Therefore, in this 
paper we use the cross-correlation betwe en the 2 1 -month Fer mi - 
LAT EGB maps with the WMAP7 maps dKomatsu et"aL[[20ICl) of 
the CMB to look for the ISW signal. 

Our analysis is also sim ilar to that of [Boughri et al.l ( ll998[) and 
[Boughn & Crittenden[ ( l2004l) since, like in their case, we also cross- 
correlate two diffuse signals. In our case it is CMB versus EGB 
whereas in their case the CMB is cross-coiTelated with hard X-ray 
background. However, for us, due to the large eiTors caused by a 
smaller number of photons, we do not expect to be able to put use- 
ful constraints on the cosmological constant or on the bias of y-ray 
sources. Instead, we will use the cross-correlation analysis to con- 
strain the nature of the EGB in the framework of the concordance 
Cold Dark Matter model with a cosmological constant (ACDM). 
Here we adopt the best-fit WMAP model (Komatsu et alj 201^ 
in which a^lT^ = 0.02267, a^h- = 0.1 131, T = 0.084, h = 0.705, 
As =2.15 X lO-** at /to = 0.05 Mpc-',and = 0.968. 

The layout of the paper is as follows: in Section|2]we briefly 
review the theoretical description of the cross-correlation analysis. 
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including the ISW signal. In Section[3]we present the various maps 
used to compute the two-point ACF and CCF. We describe the sta- 
tistical estimator used to evaluate the auto- and cross-correlation 
functions in Section|4] The results of the cross-coiTelation analysis 
are presented in Section |5] analyzed in Section|6]and discussed in 
Section|7] The conclusions are presented in Section[8] 



2 THEORY 

Our theoretical formulation of the mean y-ray emission contributed 
by unresolved sources, their auto-correlation and cross-coiTelation 
with the angular dis tribution of different type s of extragalactic ob- 
jects follows that of lAndo & Pavlidoul ( |2009|) . The treatment of the 
cross-correlati on with the CMB ma ps a nd the related IS W theory 
follows that of lBoughn et"al] ( 1 19981) and lXia et al] j2009l) . 



2.1 Mean Intensity 

The mean differential y-ray energy flux due to a population of 
sources j characterized by a y-ray luminosity function (t>j{Ly,z) 
is; 



dE 
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ATI 



J JL, 



<i>j{Ly,{l+z)E,z)LydLy 



{l+z)H{z) 



(1) 

where E is the energy, z is the redshift, H{z) = Hq[(1 +z)^Si.M + 
Q./\] is the expansion history of a flat universe with a cosmolog- 
ical constant A and the (1 term represents cosmological 

dimming of the photon energy. Ly is the source luminosity, given 
throughout the paper in erg s^', while we measure dlj/dE in 
erg cm^^s^'sr^^. Py(z) = / ^{L)LdL represents the comoving y- 
ray luminosity density of the sources at redshift z- The integration 
limits are set by requiring that the EGB consists entirely of sources 
below the photon flux detection limit S\\m (ph cm~^s~') In thi s 
case Lmax{z) = 47ld|(z)Fiin,(l + z)-2+r; jGhisellini et al]|2009l) . 
where di is the luminosity distance in the adopted cosmology and 
Tj is the photon index of the source population energy spectrum, 
assumed to be a power law (see also below). Fy^ is the energy 
flux detection limit (erg cm^^s^') , which is related to the photon 
flux 5iini through F\i^ = Sn^E, x ( 1 — ) / ( 2 — ) , where E, is the 
energy threshold of integration, typically 100 MeV. Lmin can be, 
in principle, set to zero if the total luminosity density Py(z) were 
convergent. Unluckily, for the source classes discussed below the 
extrapolation of the luminosity function to small values of Lmin 
gives a divergent total luminosity. However, an effective Lmin> if 
not given a priori, can be nonetheless computed by extrapolating 
the observed <I>j (L) and by requiring that the population of sources 
j contribute a fraction /; of the total EGB: 
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/lOO dE E Jioo dE E Jo dSm 
where dNjdS\oQ represents the differential counts of j sources and 
the factor of accounts for the units of the dNj/dSioo func- 
tion tygicanygivenjri terms of photon flux rather than energy 
flux jAbdo etaUbOlOdh . The conventional threshold of 100 MeV 
has been chosen to define the integral flux above 100 MeV 5ioo. 
This expression assumes that t he detection efficien cy can be mod- 
eled as a step function. In fact lAbdo et alj j2010dh found that the 
Fermi-LAT detection efficiency for E > 100 MeV drops as S^^ be- 
low 5iini = 3 X 10~^ ph cm~^s~'. We verified that, varying Sy,^ 
in the reasonable interval 10^^-10^^ ph cm^^s^', this stepwise 



approximation does not affect significantly the resulting redshift 
distribution of unresolved sources. From Eqs. [T] and |2] then it is 
possible to compute the mean differential flux and the fraction /j of 
the EGB emission contributed by any class of unresolved sources 
once their luminosity function and number counts are modeled the- 
oretically or extrapolated from observations. 

In this work we will assume that unresolved sources have 
power law energy spectra I{E)dE oc E^^^'dE and photon index 
Fj > 1. In this case the energy dependence drops out of the integral 
in Eq.[T]and the integrated energy flux becomes: 
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(3) 



In this paper we deal with maps of photon counts rather 
than energy flux; the photon flux (above energy E) being sim- 
ply (2 — Fj)/(1 — Fj) X Ij{> E)/E. We will considered integrated 
fluxes with three energy thresholds: /(> E = IGeV), /(> E = 
3GeV) and /(> E = 30GeV) and three possible contributors to the 
EGB: two types of blazars. Flat Spectrum Radio Quasars (FSRQs) 
and BL Lacertae (BLLacs), and star-forming galaxies. The main 
characteristics of these three populations, summarized in Table [T] 
are: 

(i) FSRQs represent a sub-class of blazars, i.e. AGNs with a 
relativistic jet pointing close to the line of sight. For the classi- 
fication of blazars as FSRQs we rely on the criteria adopted by 
lAbdo et alj | |2009|) . These authors also compute the FSRQs lumi- 
nosity function for £ > 100 MeV in three different redshift inter- 
vals, and we adopt their description under the assumption that its 
shape does not depend on the energy band. The F SRQs number 
counts have been measured by lAbdo et al. I ( l2010dl) from a larger 
sample of objects. The slope of the counts flattens in the faint 
end, implying that FSRQs at most contribute to /j ~ 25 % of the 
EGB. Their spectra are steeper than those of BLLacs, with an av- 
erage photon index of F '-^ 2.47. To model the FSRQs contribu- 
ti on to EGB we have derived an /j estimate from the logW-log5 
of lAbdoetal] ( l2010d) . We then use the luminosity function and, 
as anticipated, we enforce an effective luminosity cut Lmin- The 
calculation gives Lmin = 0.2 x lO^'^erg s^K The resulting redshift 
distribution of the flux is shown in Fig.[T] with the apparent piece- 
wise behavior resulting from the different luminosity functions in 
the different redshift bins. We stress that the value of Lmin de- 
rived in this way is not fully consist e nt wit h the minimum FSRQ 
luminosity measured by lAbdo et al] l l2009h . This affects the red- 
shift distribution in the lowest redshift bin, which however makes 
only a subdominant contribution to the total flux. Given the large 
uncertainties in the modeling of the FSRQs luminosity functions 
and their counts, this is likely a reasonable level of approximation. 
Finally, for these sources we adopt the redshift-d ependent AGN bi- 
asing function proposed by iBonoli et al.l | |2009|) in the framework 
of the semi-analytic models of AGN-black holes co-evolution: 
byiz) =0.42 + 0.04(1 +z)+0.25{l+zf. 

(ii) BLLacs are another sub-class of blazars, on average less- 
bright than FSRQs. As in the pre vious case we adop t the BLLacs 
luminosity function m easured b y 'Abdo et al.| ( |2009|) . the number 
counts determined by (^doeTal. (2010d) and enforce an effec- 
tive value of Lmin to reconcile the two predictions. The corre- 
sponding contribution to the EGB is /j ~ 12 %, while Lmin = 
6 X 10^"^ erg s~'. Finally, w e use the average photon index mea- 
sured bv lAbdo et al.l dioiOdh , i.e. F ~ 2.2 and assume that BLLacs 
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and FSRQs have the same bias factor b{z}. In this case the luminos- 
ity function is compatible with no-evolution in the different redshift 
bins and the corresponding redshift distribution in Fig. [T]has no 
piecewise behavior. 

(iii) Star-forming galaxies are fainter and much more common 
than blazars. In fact, it has been proposed that they alone could ac- 
count f or the EGB fraction which is not contributed by unresolved 
blazars. [Fields et"al] ( 1201 Ol) have shown that under the assumption 
of an Euclidean faint-end slope for the source counts, their con- 
tribution to the EGB can be as large as /j = 70 %, an assump- 
tion that we will also adopt in our analysis. Since they are very 
faint and difficult to resolve, their luminosity function cannot be 
determined experimenta lly but needs to be modeled theoretically. 
lAndo & PavlidoiJ | |2009|) have proposed a model in which the lu- 
minosity of each single source scales with the star-formation rate 
and the gas mass fraction. Since in this model the y-ray emissivity is 
rescaled from that of the Milky Way, the underlying assumption is 
that most of the y-ray photons are emitted from de-evolved versions 
of our own Galaxy. Current theoretical uncertainties and weak ob- 
servational constraints do not allow to discrimina te among the sim- 
ple model proposed bv lAndo & Pavlidoul ^2m^ and the more re- 
cent (and more sophisticated) one presented bv lFields et al.l j20IOl) . 
For this reason in this paper we have decided to adopt the first one. 
The energy spectrum of star-forming galaxies is characterized by a 
strong pionic peak at £ ~ 0.2 GeV, a feature that also determines 
the energy dependence of their contribution to the EGB. In the en- 
ergy range we are interested in (£ > 1 GeV), their energy spectrum 
is fairly well appro ximated by a power law behavior with F '-^ 2.475 
jFields et alJbOld) that allows us to use Eq.[3l Finally, we assume 
that star-formin g galaxies are unbias ed, i.e. = I, as suggested 
by observations jAfshordi et al.ll2004l) . The results of the auto- and 
cross-correlation studies do not chang e significantly as long as the 
bias of these galaxies is close to unity jAndo & Pavlidoull2009l) . 

In Fig. [T] we show the redshift dependence of the normalized 
energy flux per unit redshift t/ln/(> E)/dz for the three proposed 
EGB sources: FSRQs (red, dashed), BLLacs (black, continuous) 
and star-forming galaxies (blue, dot-dashed). This function, pro- 
portional to the integrand in Eq.|3] represents the contribution to the 
EGB from the sources in a specific redshift range. In the BLLacs 
scenario, the contribution to the EGB signal is relatively local, i.e. 
produced by a population of faint, nearby sources. In the FSRQs 
case the signal mostly comes from 2 and drops to zero at 
z < 0.9. This reflects the fact that FSRQs are rare bright objects 
which, for z < 0.9, would have L„„„ < Lmin in Eq[T] i.e. would be 
above detection threshold and removed from the map. As a con- 
sequence, the cross correlation of the FSRQs signal with catalogs 
of objects whose number density peaks at low redhifts is expected 
to be zero, as we shall se in Section|5] Star-forming galaxies rep- 
resent an intermediate case in which the signal is produced over a 
relatively broad redshift range around z ~ 1 . 
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Figure 1. Normalized y-ray flux per unit redshift t/ln/(> E) /dz as a func- 
tion of z for three different source classes: FSRQs (red, dashed), BLLacs 
(black, continuous) and star-forming galaxies (blue, dot-dashed). 



Table 1. Investigated EGB contributors. All values refer to £ > 3 GeV and 
S < Siii„. The luminosity f unction, ^ and sp ectr al pai'ameter s for FS RQs 
and BLLacs are taken from lAbdo et^ j20o4 and lAbdo et afl j2010dh . 
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We also assume that the y-ray sources trace the underlying fluctua- 
tions in the mass density according to some linear biasing prescrip- 
tion that may depend on the redshift: 



5„^(z,x) = by{z)5m{z,x) = by{z) 



Pm{z,x) - p„{z) 



(5) 



where p,„ indicates the mass density and by{z) is called the biasing 
function. 

Putting all together, the expected fluctuation in y-ray energy 
flux is: 



, l{n)-I Jil+z)-^H{z)-'pj{zMz)8,n{z,x)dz 



f{l+z)-^H{z)-'py{z)dz 



where I = I(> E) indicates the y-ray mean flux and /(n) = /(> 
£,11) is the energy flux along the generic direction n. 



2.2 Fluctuations in the y-ray flux 

To compute the predicted auto- and cross-correlation signals we 
need to model the fluctuations of the y-ray flux. These fluctuations 
arise from local deviations from the y-ray luminosity density Py(z) 
that we assume to be proportional to the deviations from the mean 
number density of sources ny{z) = / ^{L)dL: 

^ Pliz,x)-Py{z) ny{z,x) ~ ny{z) 5. , . ,., 

^^^'^'"^^ — m — ^^M^) — " ^ ^ 



2.3 Two-point Angular Correlation Function and Spectrum 

We can now compute the expressions for the two-point ACF of 
the EGB fluctuation field E and the two-point CCF between E and 
another fluctuation field obtained from a source catalog, 5. 

The general expression for the two-point angular correlation 

is 

(6/(ni)5y(n2)) = ^ ^c'/p,[cosm , (V) 
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where / and J are the two fields and the angular spectrum is given 
by: 

C'/ = \jk'P{k)[G\{k)][G'j{k)]dk. (8) 

and P{k) is the present-day power spectrum. 

For example, for the EGB ACF one has I = J = E with 



_ f{l+z)-^H{z)-'pyiz)by{z)D{z)j,[kxiz)]dz 
f{l+z)-^H{z)-^pj{z)dz 



Gf{k) 



(9) 



where ji[kx{z)] are spherical Bessel functions, D{z) is the linear 
growth factor of density fluctuations and x{z) is the comoving dis- 
tance to redshift z- In our analysis, we use the public code CAMB 
jLewis & Bridi3l2002h to generate the linear power spec trum of 
density fluctuations and the Hal of it jSmith et alj|2003h built-in 
routine for non-linear correction to obtain the fully-evolved, non- 
linear matter power spectrum P{k) at any epoch. We note that non- 
linear corrections do not affect our results significantly. 

In case of cross-correlation with a fluctuation field of discrete, 
unresolved sources one has I = E, J = S and 



Kk)-I 



dN{z 



-bs{z)D{z)ji[kx{z)]dz., 



(10) 



where dN{z)/dz and by{z) represents the redshift distribution and 
the bias factor of the sources that do not necessarily coincide with 
the y-ray sources. 

Finally, if one cross coiTelates EGB with another diffuse sig- 
nal, like the temperature fluctuation field obtained from the CMB 
maps (I = E, J = T) then 



Gj{k) = -2j 



-^hW)\dz ., 
dz 



(11) 



where <I> represents the gravitational potential. 

In this case the cross-coiTelation signal in Eq.|7]represents the 
ISW effect, expected if the EGB were contributed by sources that 
trace the underlying mass distribution. 

In Section[5]we will use Eqs.|7l[8l|9l[iO]and[TT]to predict the 
expected auto- and cross-correlation signal that will be compared 
with data. 



3 MAPS 

In this section we describe the various maps (residual isotropic 
Fermi-LAT maps, WMAP7 CMB, SDSS QSOs, SDSS, LRGs, 
2MASS galaxies, NVSS radio galaxies) that will be used for the 
auto- and cross-correlation analysis. 



3.1 Fermi-LAT maps 

The Large Area Telescope (LAT) is the primary instrument on- 
board the Fermi Gamma- ray Space Telescope launched in June 
2008 jAtwood et af]|2009l) . It is a pair-conversion telescope cov- 
ering the energy range between 20 MeV and 300 GeV. Due to 
its excellent angular resolution (~ 0.1° above 10 GeV), very large 
field of view (~ 2.4 sr) and efficient rejection of background from 
charged particles, it represents a key experiment for y-ray astron- 
omy. Fermi-LAT is continuously scanning the sky in survey mode, 
providing a complete image of the sky every 3 hours. 

The data gathered by the telescope are available online Q. For 



http://fermi.gsfc.nasa.gov/ssc/data/access/ 



the present analysis we used, unless otherwise specified, 21 months 
of data. We used only events classified as class 4. Class 4 events 
form the dataset with the lowest level of CR background contam- 
ination currently available for LAT data analysis. Deta ils on this 



event classification are described in 



Abdo et al 



(l2010d) . We also 

use events labeled as class 3 (cf. lAtwood et all ( |2009|) )" for consis- 
tency checks0. Class 3 events have a larger residual contamination. 
We apply a cut of 100° on the zenith angle, 52° on the satellite 
rocking angle and 65° on the inclination angle to reduce the con- 
tamination from Earth al bedo. The counts t hen have been pixelized 
in the HEALPix format jGorski et al.l2005h with Njide = 64, corre- 
sponding to Wpix =49152 pixels with an angular size 0.92° x 0.92°. 
The majority of the se events come from the resolved point sources 
jAbdo et allboiObl) and from the Galactic diffuse emission due to 
y-rays produced in the interaction of cosmic rays with the interstel- 
lar medium. Both need to be removed to extract the EGB signal. 

To account for the Galactic diffuse foreground we adopted the 
model gll_iem_v02 . fit0. In the following this will be re- 
ferred to as VI. The model is based on fits of the LAT data to tem- 
plates of the HI and CO gas distribution in the Galaxy as well as an 
Inverse Compton model obtained with th e GALPROP code an d 
a further template for the Loop-I region (ICasandiian et al. lliooi). 
The model, based on 1 year of data taking by Fermi-LAT, describes 
well the Galactic diffuse emission over the sky, with some excep- 
tion, most prominently in th e regions whic h have been associated 
with giant gamma-ray lobes jSu et al.l2010l) . We will thus also use 
a preliminaiy version of a refined model in development, which, 
further, includes templates specifically for these structures, sensi- 
bly reducing the residuals in these regions. This is model V2. In 
the next section we check the impact of these models on the ACFs- 
CCFs and compare the results with those obtained by masking out 
the areas where the Bubble and Loop-I emission are more promi- 
nent. The mask is obtained from the tabulated contours given in 
dSu et alj|2O10l) excluding the whole region enclosed by Loop-I in 
the North Galactic sky. In addition, for both models, we masked out 
a 20° strip above and below the Galactic plane where the EGB sig- 
nal is largely subdominant. And in order to check the robustness of 
our results, we vary the strip mask width from b = 20° lo b = 50°, 
in steps of Afc = 10°. 

Point so urces in the 1 ye ar Fermi-LAT point-source cata- 
logue IFGL jAbdo et alj|2010bl) are also masked as described in 
ICuoco et al] ilOld) . The point source masking is adaptive, in the 
sense that more intense sources, more likely to bias the results, are 
masked with a larger circle (up to 2-3 degree radius), while fainter 
ones are masked with a smaller circle down to about 1 degree (see 
Fig.llJ. 

To subtract the Galactic foreground from the data we have 
first generated the energy dependent exposure maps with the 
gtexpcube routine available in the public Fermi-LAT science 
toolsQ Version p6v3 class 4 of the LAT instrument response func- 
tions (IRFs) and the same cuts used for the event selections were 
used for the exposure map generation. Maps of expected fore- 
ground counts from our two adopted Galactic diffuse emission 



^ Class 4 events are also referred to as "dataclean" class while the 
union of class 3 and class 4 events is also referred to as "diffuse" class. 
See |http://femu. gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone I 
/Cicerone_Data_ExpIoration/Data_preparation.html 

' A more detailed description can be found at 
|http://fermi. gsfc.nasa.gov/ssc/d ata/access/lat/BackgroundModels. html I 
^ http://galprop.stanford.edu/ 
' See http:// fermi.gsfc.nasa.gov/ssc/ 1 
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Figure 2. Residual gamma-ray emission obtained from analyzing 21 
months of Fermi-LAT data for energies E > I GeV (top), E >3 GeV (sec- 
ond from the top), £ > 30 GeV (second from the bottom) and the mask cov- 
ering the point sources in the 1 year Fermi-LAT point sources catalog and 
the region with \b\ < 20° (bottom panel). The two contours in the mask in- 
dicate the regions which are further removed when the Lobes/Loop-I mask 
is applied. The residual maps have been smoothed with a 2° Gaussian filter 
to remove small scale Poisson noise. All maps are in Galactic coordinates 
and have a resolution Wsidc = 64 (~ 0.92°). All Fermi-LAT maps have been 
cleaned by removing the V2 Galactic model and all multipoles with £ 10 
(see text). 



models have been calculated by applying the exposure maps to the 
models and convolving them with the Point Spread Function (PSF) 
of the LAT averaged over the fie ld of view and the relev ant energy 
range. The GaRDiAn package jAckermann et al.ll2008[) was used 
for this task. We thus include explicitly the PSF convolution, al- 
though we stress that in the energy range considered the PSF width 
is always smaller or at most comparable to our minimum angular 
bin (approximately 1°). However, further checks of the effect of the 
PSF on the ACF are considered in the next section. After the sub- 
traction, the residuals should be dominated by the EGB signal and 
by the residual isotropic instrumental background. We have pro- 
duced three maps containing events with energy E > I GeV, E > 3 
GeV and £ > 30 GeV, respectively. The energy cut was performed 
to reduce the chance of residual contamination from Galactic emis- 
sion and to evaluate different trade-offs between clean maps and 
reasonably large statistics. 

Since models VI and V2 are not peifect, residuals can contain 
some spurious signal on large scales that may affect the correlation 
studies. The impact on the CCF is small since these spurious resid- 
uals are not expected to correlate with the extragalactic signal, but 
might be significant on the ACF. For all the maps, we thus always 
remove the residual monopole and dipole contribution. This is per- 
formed using the specific HEALPix routine, remove_dipole. 
In addition, we have implemented a more aggressive cleaning pro- 
cedure in which, besides removing the monopole and dipole, we 
have expanded the residual maps into spherical harmonics using 
the HEALPix tool anaf ast and have removed all the multipoles 
up to ^ = 10 that contribute to the large scale signal. We will refer to 
the residual maps with all £ < 10 multipoles removed as £10 maps 
whereas £1 indicates the maps in which only dipole an monopole 
have been removed. The validation of this procedure and testing 
for possible, undesired systematic effects is described in the next 
section. 

The three £10 residual counts maps are shown in Fig. [alto- 
gether with our fiducial mask removing the region within 20 de- 
grees from the Galactic Plane and the point sources. To better il- 
lustrate the fluctuation properties of these maps, the mean has been 
removed, thus showing only fluctuations around zero. Furthermore, 
the maps have been smoothed with a 2° Gaussian filter to remove 
small scale Poisson noise. 



3.2 Validation and checks 

The scope of this section is to describe the various tests that we 
have performed to check and validate our cleaning procedure, to 
reveal the presence of possible systematic effects, to correct them 
and estimate their impact on the ACF and CCF analysis. For brevity 
we only show a subset of the results obtained, but we stress that all 
robustness tests described in this section have been repeated for 
each of the ACF/CCFs analyses presented in this paper. 

3.2.1 Validating the £10 cleaning procedure 

To validate the £10 cleaning procedure we have applied this clean- 
ing technique to a set of mock data mimicking the characteristics 
of the real 21-month Fermi-LAT data. The mock, Monte Carlo- 
generated data-sets are simulated with the gtobssim routine us- 
ing the p6_v3_dif f use Instrument Response Functions and in- 
cluding the contribution from three signals (i) a Galactic diffuse 
component generated using the VI Galactic model, (ii) an extra- 
galactic diffuse and isotropic component (hence with no intrinsic 
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CCF with the CMB (right) relative to V2 model and £ > 3 GeV cut. 
Top panels show the impact of the different cleaning procedures for 
the V 1 galactic model. Middle panels are relative to the V2 galactic 
model. Bottom panels show the impact of different \b\ cuts. 

The results of these tests can be summarized as the following: 
(i) for the ACT analysis the Galactic model V2 performs signif- 
icantly better than VI and therefore will be used throughout this 
work, (ii) The ACF is quite sensitive to different cleaning pro- 
cedures and both the Bubbles/Loop-I mask and the HQ cleaning 
are required to get consistent results. After this cleaning, the auto- 
coiTelation signal is stable for Galactic cuts \b\^ 20° . For these rea- 
sons, we will consider (and show) only the ACFs computed using 
the V2 model with f 10 cleaning, Bubbles/Loop-I mask and a Galac- 
tic cut \b\ = 20°. (Hi) The CCFs are remarkably insensitive to clean- 
ing procedures and masks already when using the V 1 model, apart 
the case of no Bubbles/Loop-I mask and no ^10 cleaning which, not 
including any modeling of the Bubbles, leaves significant residuals 
which seem to bias the CCF result. Overall, this test show that the 
ACF and CCF signals converge to a stable result when progres- 
sively more aggressive cleaning procedures and latitude cuts are 
applied. For consistency with the ACF case we will consider (and 
show) the CCFs computed using the same Fermi-LAT maps used 
for the ACF analysis, although we did check for the robustness of 
all the CCFs with respect to the different cleaning procedures and 
masks. 



Figure 3. ACFs (left panels) and CCFs with the CMB (right panels) for dif- 
ferent cleaning procedures {£10 and £1, with/without Bubbles/Loop-I masft) 
for the VI Galactic model (top panels), V2 Galactic model (middle panels), 
and for different Galactic latitude cuts (bottom panels). All plots refer to 
Fermi E > 3 GeV residual maps. 

ACF) generated with the isotropic_iem_v02 . txt0 model 
and ( Hi) a signal contributed by a population of AGNs generated by 
sampling their observed logA'-logS distribution ( I Abdo et al 
of which we mask the 1200 brightest ones. No intrinsic clustering 
is assumed for the simulated AGNs which are distributed across 
the sky randomly. The same pipeline described in the previous 
section is then applied to the mock datasets to obtain the residual 
maps from which ACFs and CCFs with the CMB are calculated. 
We observe that the resulting ACFs and CCFs are generally nicely 
consistent with zero. Some systematic features at the level of la 
however appear as, mainly, a negative offset of the ACF for the £10 
maps. Indeed, we also observe this offset for the ACF of real data. 
We thus plot a gray region around the measured ACF to indicate 
the presence of this systematic. 

3.2.2 Robustness to the cleaning procedure, Galactic model and 
masking 

To test the robustness of our results with respect to the cleaning 
procedure {£1 vs. ^10), Galactic model (VI vs. V2) and Galactic 
mask {\b\ < 20°, 30°, 40°, 50° and Bubbles/Loop-I cut) we have 
computed the residuals of real Fermi maps using the different com- 
binations of these procedures at all energy cuts. Then we have com- 
puted and compared the various ACFs and CCFs with all catalogs 
described in the following section. An example of these tests is 
shown in Fig. [3] The different panels show the ACF (left) and the 

* http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModeIs.htmI 



3.2.3 Robustness to the event class type 

A further check is performed comparing the ACF and CCFs pre- 
viously obtained from class 4 events maps with those computed 
from the class 3 and class 4 events together. In this second case the 
increased contamination from residual isotropic instrumental back- 
ground is not expected to affect the CCFs although, of course, it 
can make the error bars slightly larger, which is indeed what was 
found. For the ACFs more effects come into play. For example, the 
enhanced isotropic instrumental component is expected to decrease 
the amplitude of the auto-correlation signal. However, since the 
measured ACF is consistent with zero, the only effect is, again, to 
slightly enlarge the error bars, which is indeed observed. Of course, 
the additional instrumental background associated to class 3 events 
may not be completely isotropic. However, the fact that the mea- 
sured ACF is consistent with that of the class 4 events indicates 
that possible deviations from isotropy are too small to affect our 
analysis. 

3.2.4 Sensitivity to the PSF of the detector 

The angular resolution of the detector, characterized by its PSF, 
in principle introduces a correlation between neighboring pixels 
which may increase the auto-correlation signal on angular scales 
comparable with the PSF of the detector (causing the so-called am- 
plification bias). The effect, however, is expected to be small since 
our bin-size is significantly larger than the original resolution of 
the maps. To estimate its impact on our analysis we have first com- 
puted the window function in multipole space associated to the PSF, 
given by W,{E) = f d8P/(cos(e))PSF(E,e), where P/(cos(e)) are 
the Legendre Polynomials and PSF(E,0) is the energy dependent 
PSF. We have then multiplied W/ for the expected, intrinsic cor- 
relation given by Eq[8]and integrated back to configuration space. 
The result is shown in Fig.|4]in which we compare the expected in- 
trinsic ACF (continuous, black curve) to that convolved with the 
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Figure 5. The 7-yeai' WMAP ILC map in Galactic coordinates with the 
resolution A'sye =512. 
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Figure 4. The effect of finite PSF on the ACF of BLLacs. 



PSF at £• =1 GeV, 3 GeV, and 30 GeV which has a 68% con- 
tainment radius of 0.8° (green, dot-dashed) ,0.3° (blue, dotted), 
0.1° (red, dashed) respectively ( I At wood et al. I I2OO9I) . Notice also 
that we compared the results with a simple Gaussian approxima- 
tion for the window function, = [exp(— where o^, is 
the width of the beam, finding almost indistinguishable results. As 
can be seen, only for £ ^ 1 GeV some small effects can be appre- 
ciated, while for E >3 GeV this amplification bias is completely 
negligible. 
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3.2.5 Robustness to the event conversion type 

As a further robustness te st we also compute d the ACF of those 
events labeled as front jAtwood et al.ll20()9l) which are photons 
converting in the front part of the detector and have a significantly 
better PSF with respect to the rest of the events which, instead, 
convert in the back part of the detector where thicker converter foils 
increase the chance of large-angle scattering which deteriorates the 
tracking accuracy. They amount to about half of the events detected 
by the Fermi-LAT detector. We performed this test to check the 
robustness of the CCF with SDSS-LRGs and 2MASS galaxies in 
Sections 15.51 and \5l6\ i.e. in the two cases in which we find some 
features in the CCFs. Apart from some increase in the error bars 
due to the halved statistics, the results did not change significantly. 



3.3 WMAP7 ILC 

In order to search for the ISW effect, we cross-correlate the Fermi- 
LAT 21 -month EGB maps with the CMB maps derived from the 
7-year WMAP data. More precisely, we use the WMAP Internal 
Linear Combi nation (ILC) map w ith N^ide = 512 provided by the 
WMAP team jKomatsu et alj201Cl ). shown in Fig.[31 This ILC map 
was already built to minimize the Galactic and other foreground 
contaminations. F or the WMAP map, we use the "KQ75" mask 
jGold et al.ll2oToh corresponding roughly to the "KpO" cut in the 
3-year data release. In our calculations, we downgrade it to the res- 
olution A'side = 64, to match that of the Fermi-LAT maps, and set 
the weight wx = for all pixels including at least one m asked high 
resolution pixel jRaccanelh et al.ll2008|jxia et alj|2009h . 



Figure 6. Normalized redshift distributions, dN/dz, of the different types 
of objects considered for our cross-conelation analysis. SDSS DR6 QSOs 
(black, continuous Une), 2MASS galaxies (red, dashed), NVSS galaxies 
(blue, dot-dashed) and SDSS DR6 LRGs (cyan, short-dashed). 



3.4 Discrete sources maps 

One of our goals is to cross-correlate the EGB maps with different 
classes of sources that trace, but not necessarily coincide with, the 
EGB sources. Since all luminous objects trace, with a different de- 
gree of bias, the same underlying distribution of matter, it makes 
sense to cross correlate the EGB with the following sources; Op- 
tically selected quasars, ii) luminous radio galaxies, «') IR-selected 
galaxies, and iv) LRGs. (;) and (/;) span the same, broad redshift 
range as the FSQRs whereas (Hi) and (iv) span much narrower red- 
shift ranges that overlap with those of BLLacs and Starforming 
galaxies. Below, we provide some details on the different catalogs 
considered in our analysis. In Fig|6]we show the redshift distribu- 
tion, dN/dz, of the four catalogues we have considered in this pa- 
per. All distributions are normalized to unity. 2MASS galaxies and 
LRGs trace the large scale structure of the local universe and, from 
Fig.[T] we see that we can expect some cross-correlation signal only 
if the bulk of the EGB is contributed by star-forming galaxies or 
BLLacs. On the contrary, in the case of QSOs and NVSS galaxies, 
a positive cross-coiTelation may be expected if EGB were prefer- 
entially contributed by a population of high-redshift objects like 
FSRQs. However, the broad redshift distribution of these objects 
might also allow to pick up some cross-correlation signal provided 
by a population of low-redshift y-ray sources. 
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Figure 7. The number count map of SDSS DR6 quasar catalogue in Galac- 
tic coordinates with the resolution A'sidc = 64. 



Figure 8. The number count map of 2MASS extended source catalogue in 
Galactic coordinates with the resolution A'side = 64. 




3.4.1 SDSSDR6QS0 

We us e the SDSS DR6 quasar catalog released by [Richards et alj 
( I2OO9I) (hereafter DR6-QSO). This catalog contains about N^^^lif 
quasars with photometric redshifts between 0.065 and 6.075, cov- 
ering almost all of the northern hemisphere of the Galaxy plus three 
nan'ow stripes in the southern, for a total area of 8417deg2 (~ 20% 
of the area of the whole sky). The DR6-QSO data se t extends previ- 
ous s i milar SDSS data s ets with ~ 95% efficiency jRichards et al.l 
l2004lMversetal.ll2006l) . The main differences are due to the fact 
that DR6-QSO probes QSOs at higher redshift and also contains 
putative QSOs flagge d as to have ultra violet excess (UVX objects). 
We refer the reader to lRichards et al.l(l2009l) for a very detailed de- 
scription of the object selection with the non-parametric Bayesian 
classification kernel density estimator (NBC-KDE) algorithm. 

We rely on the electronically published table that contains 
only objects with the "good" flag with values within the range [0,6]. 
The higher the value, the m ore probable for the object to be a real 
QSO (Ri chards et al.ll2009l) . We only consider the quasar candi- 
dates selected via the UV-excess-only criteria "uvxts=l", i.e. ob- 
jects clearly showing a UV excess which should be a signature of a 
QSO spectrum. We are left with A'q ~ 6 x 10^ quasars. In Fig.|7]we 
show the number counts map of the SDSS DR6 quasar catalogue 
in Galactic coordinates. 

In order to determine the mask of the actual sky coverage 
of the DR6 survey, we generate a random sample with a suffi- 
ciently large number of galaxies using the DR6 database to ensure 
rou ghly uniform sa mpling on the SDSS CasJobs website. Follow- 
ing Kia_et_al, liooi), besides the pixel geometry mask, we also add 
the foreground mask by cutting the pixels with the g-band Galactic 
extinction Ag = 3.793 x £(6 — K) > 0. 18 to account for reddening 
that is the main systematic effect. 

The redshift distribution function dN jdz of the DR6-QSO 
sample is approximated by the function: 



dN 
dz 



P 



rf m+l \ _m+l 



exp 



^0 



(12) 



where three free parameters are m = 2.00, P = 2.20, and zq = 1.62 
jXia et al.ll2009l). W e choo s e a constant bias = 2.3 as found by 
iGiannantonio et al.l l l2008h : IXia et al.l l l2009h to calculate the theo- 
retical prediction from the best-fit WMAP model adopted in this 
work. 



3.4.2 2MASS 

We use the 2 Micron All-Sky Survey (2MASS) extended source 
catalogue jjarrett et all l2000r) . 




Figure 9. The number count map of NVSS radio sources in Galactic coor- 
dinates with the resolution A'side = 64. 



ies with mean redshift (z) ~ 0.072, as shown in Fig. [8] We se- 
lect galaxies according to their A^s-band isophotal magnitude K2Q, 
measured inside a circular isophote with surface brightness of 20 
mag/arcsec^. These magnitudes are corrected for Galactic extinc- 
tion using the infrared reddening maps: ATjq = A'20 — Ai^, where 
the extinction Aj^ = 0.367 x {B ~ V) . In our analysis, we use the 
flux cut 12.0 < < 14.0. We only include objects with a uni- 
form detection threshold (use_src = 1), and remove known arti- 
facts (cc_flag 7^ a and cc_flag 7^ z). Furthermore, we exclude areas 
of the sky with high red dening using the infrared reddening maps 
bv lSchlegel et al.l jl998l) . discarding pixels with A^ > 0.05, which 
leaves approximately 67% of the sky unmasked. 

In this case, the free parameters of the re dshift distribution in 
Eq. | 1 2lare m = 1.90, P = 1.75, and zo = 0.07 jGiannantonio et aT] 
200 8l) . while as constant bias we use bs = 1.4 as found by 



Rassat et all j2007l) . 



3.4.3 NVSS 



which contains ~ 770000 galax- 



The NRAO VLA Sky Survey (NVSS) jCondon et alll 19981) offers 
the most extensive sky coverage (82% of the sky to a completeness 
limit of about 3 mJy at 1.4 GHz) and contains 1.8 x lO*" sources. 
Here, we include in our analysis only NVSS sources brighter than 
10 mJy, since the surface density distribution of fainter sources suf - 
fers from declination-dependent fluctuations jBlake & Wallll20()2h . 
We also exclude the strip at \b\ < 5°, where the catalog may be sub- 
stantially affected by Galactic emissions. The NVSS source surface 
density at this threshold is 16.9 deg^^. 

The reds hift distribution at thi s flux limit has been recently 
determined bv lBrookes et al.l ( |2008|) . Their sample, complete to a 
flux density of 7.2 mJy, comprises 110 sources with 5 ^ 10 mJy, 
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Figure 10. The number count map of MegaZ LRG catalogue in Galactic 
coordinates with the resolution NsHs = 64. 

of which 78 (71%) have spectroscopic redshifts, 23 have redshift 
estimates via the K — z relation for radio sources, and 9 were not 
detected in the A^-band and therefore have a lower limit to z- We 
ad opt here the smooth d escription of this redshift distribution given 
bv ldeZottietai] ( l201(]|) : 

^(z) = 1.29 + 32.37z-32.89z2 + 11. 13z^- 1.252"* . (13) 
dz 

Here, we simply use a constant bias Z^s = 1.5 to calculate the theo- 
retical prediction. 

3.4.4 SDSS DR6 LRG 

We use the updated MegaZ LRG sample jAbdalla et al] |2008|) , 
which contains ~ 1.5 x 10*' galaxies from the SDSS DR6 in the 
redshift range 0.4 < z < 0.7 with limiting magnitude i < 20, as 
shown in Fig. (TO) To reduce stellar contamination, there is a vari- 
able of the MegaZ neural network estimator 5sg, defined such that 
5sg = 1 if the object is a galaxy, and Sjg = if it is a star. For a con- 
servative analysis, we choose a cut Sjg > 0.2, which is reported to 
reduce stellar contamination below 2% while keeping 99.9% of the 
galaxies. In addition to the SDSS DR6 geometry mask, we also add 
two foreground masks to account for seeing (removing pixels with 
median seeing in the red band larger than 1.4 arcsec) and reddening 
(removing pixels with m edian extinction in the red band Ai- > 0. 1 8) 
jGiannantonio et alj20 0S) . 

The redshift distribution function in this case is found di- 
rectly from the photometri c redshifts that are given in the cata- 
logue (Collisteret al. 2007t), while we set the constant bias to 
1.8 tOiannantonio et al..2008l) when calculating the theoretical pre- 
dictions. 



4 CORRELATION ANALYSIS 

In order to calculate the ACF and the CCF all maps have been re- 
binned to a resolution of N^ide = 64. To estimate the ACF of the 
residual Fermi-LAT maps we use the following estimator c*(9): 

£''i^) = J^^Li"i-")i"J-''i)' (14) 

where n[ is the number of y-ray photons in each pixel and ti is the 
mean photon counts corresponding to the average EGB signal. The 
sum runs over all the pixels with a given angular separation. For 
each angular bin centered around 6, A'e is the number of pixels 
pairs separated by an angle within that bin. The CMB maps consist 
of temperature difference maps. Therefore we replace the n; and ii 



with the temperature T; in each pixel and the average temperature 
f of the CMB map. The CCF estimator (in /jK) between the Fermi- 
LAT map and the CMB temperature map c*^(6) reads: 

e"{e) = r^^U'r<-T){nj-n). (15) 

Moreover, we are also interested in the cross-correlation between 
Fermi-LAT maps and the different maps of discrete sources. In this 
case, the CCF estimator is c's(9): 

£'m = ^L('^i~"r-^^- (16) 

Since we are using a resolution of A'side = 64, for which the pixel 
size is 0.92°, in our calculation we use A'bin = 5 linearly spaced 
angular bins in the range 1° ^5 ^ 9°. 

To compute the errors, we estimate the covariance matrix of 
the ACF (CCF) data points using jackknife resampling method 
jScranton et al .I2OO2I) . This method divides the data into M patches, 
then M subsamples are created by neglecting each patch in turn. 
These patches have roughly equal area. In practice, we firstly list 
the whole set of pixels covered by the survey, and then divide them 
into M = 30 patches that do not have very similar shape, but have 
roughly equal area (i.e. equal number of pixels). The covariance 
estimator is: 

M- 1 

Cij = -JT I {cm-C{Qi)\ [Q(9y) -C(e,0] , (17) 

A=l 

where Ca.(6,) are the observed ACF (CCF) of the M subsamples 
in the ;-th angular bin and C(6i) are the mean ACF (CCF) over M 
realizations. The diagonal part of these matrices gives the variance 
of the ACF (CCF) in each bin C\ = af, while the off-diagonal part 
represents the covariance between the angular bins. We also change 
the number of patches M and verify that the covariance matrix is 
stable. 



5 RESULTS 

In this section we compute the ACF and the CCF of the different 
map combinations and compare the results with model predictions. 
The ACF signal could constrain the nature of the sources that con- 
tribute to the EGB. The presence of the ISW effect would quantify 
how well the EGB sources trace the underlying mass distribution 
and reveal the presence of a cosmological constant term. Finally, 
the strength of the CCF signal would indicate how closely a given 
class of objects trace the sources of the EGB. For each correla- 
tion analysis we will show measured quantities and theoretical pre- 
dictions. The comparison between the two is only performed in 
configuration space, i.e. we over-plot the theoretical ACF and CCF 
(continuous curves) with the data points with error bars. In addi- 
tion, we also will show mode predictions in Legendre space, i.e. 
through the expected angular power spectra, since this is the direct 
prediction of the model, and model to model comparison is easier 
in Legendre space. All theoretical predictions use Eq|8]to model 
the correlation between the y-ray signal supposedly contributed by 
a single type of objects (FSRQs, BLLacs or Star forming galaxies) 
and a second fluctuation fields (the y-ray background for the ACF, 
the CMB for the ISW signal and the discrete source catalogs for the 
cross-coiTelation analyses). 
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5.1 Auto-correlation analysis 

In Fig[TT]we show the ACF of the residual maps obtained using the 
V2 Galactic model with l\Q cleaning, Bubbles/Loop-I masking and 
a Galactic cut \b\= 20° . Different symbols indicate different energy 
cuts: black square is for E > \ GeV, red dot is for £ > 3 GeV and 
magenta triangle is for £ > 30 GeV. Results for the £ > 30 GeV 
case are shown in separate panel since the error bars are signifi- 
cantly larger, lo error bars were computed using the jackknife pro- 
cedure. The three continuous curves in this plot and in all plots of 
the cross-correlation functions discussed in the following sections 
represent the theoretical predictions from Section l23] obt ained as- 
suming that each source class (BLLacs, FSRQs or Star Forming 
galaxies), specified by their luminosity density distribution Py(;), 
contribute to a fraction of the EGB, fj, listed in Table[T] 

At small angular separations (9 < 2°) the auto-coiTelation sig- 
nal is consistent with zero with all energy cuts but £ > 1 GeV. 
Theoretical models do predict a weak auto-correlation signal at 
these angular separations that can be regarded as the typical an- 
gular size of the y-ray emitting element. However, the predicted 
auto-correlation is much weaker than the measured one, at a level 
that would be indistinguishable from zero with the current uncer- 
tainties. Although a correlation signal at this level would still be 
possible from a contribution of unresolved point sources, we find 
that the signal is not very robust to the different cleaning methods 
and progressively disappears when we apply larger \b\ cuts. We 
therefore regard this as a spurious feature of the ACF at low ener- 
gies where the Galactic contribution is stronger and its subtraction 
more prone to systematic errors. At larger separations (6 > 4°) the 
ACF is slightly negative. This is a spurious features induced by the 
/lO cleaning procedure as we have shown the Monte Carlo analysis 
presented in Section [3.2. 11 The corrected signal would be consis- 
tent with zero. From the Montecarlo we estimate that this system- 
atic offset is at the level of the lo statistical eiTor. We thus show in 
the plot also a systematic uncertainty band obtained doubling the 
lo statistical uncertainties. 

To highlight the differences among the different source classes 
we plot the angular auto-power spectra in Fig. [12] BLLacs and star- 
forming galaxies have similar spectra, with more power on large 
scales than the FSRQ model. The situation is reversed at small 
scales. This difference simply reflects the fact that in the first two 
cases the y-ray emission peaks at moderate redshifts z = [0.5, 1] 
while for FSRQs the bulk of the y-ray signal is produced at z > 2, as 
shown in Fig.[T] Since all models trace the same mass density field, 
i.e. assume the same P{k), the power shift in Fig.lll [reflects the fact 
that the same physical scale is preferentially seen at different an- 
gles in the different models: large angles (small () for BLLacs and 
star-forming galaxies that typically sample the universe at low red- 
shifts; small angles (large £) for FSRQs that preferentially samples 
the universe at high redshifts. The larger amplitude of the spectrum 
for BLLacs to the spectrum of star-forming galaxies simply reflects 
the different bias factors of the two classes of objects. 



5.2 ISW 

In Fig.[T3]we show the CCF of the EGB with the CMB, i.e. the ISW 
signal. The symbols refer to the CCF estimated from the 21 -month 
Fermi-LAT maps and the WMP7 ILC map and the continuous lines 
represent model predictions. Different symbols indicate different 
energy cuts. 

We note that the expected CCF signal is positive out to large 
angular separations and is only a factor 3-4 smaller than the width 
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Figure 11. ACFs estimated from the 21 -month Fermi-LAT EGB map for 
\b\ > 20° in three energy bands: £ > 1 GeV, i? > 3 GeV (upper panel) and 
£ > 30 GeV (below panel). Model predictions for different types of sources 
are represented by continuous curves: FSRQs (black, continuous), BLLacs 
(red, dashed) star-forming galaxies (blue, dot-dashed). The gray band indi- 
cates the systematic uncertainty coming from the foreground cleaning pro- 
cedure estimated to be approximately equal to the 1 a statistical errors (see 
text). In the top panel only the band for the 1 GeV case is shown. 
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Figure 12. Angular auto-power spectra for different EGB source models . 
Different line-styles characterize different models. Black, continuous curve: 
FSRQs; red, dashed: BLLacs; blue, dot-dashed: star-fonmng galaxies. 

of the error bars. This has to be compared with the ACF case where 
the expected signal is, instead, more than a factor of 10 smaller 
than the the experimental error bars. This result, which reflects the 
better statistics of the CMB maps, indicates that the goal of detect- 
ing the ISW signature in the EGB is not unrealistic. However, with 
the limited statistics of the 21 -month data the ISW signal is con- 
sistent with zero with all energy cuts. This (null) result is robust to 
different cleaning methods and latitude cuts. 

From a theoretical point of view we notice that the ISW sig- 
nal is expected to be larger in the BLLacs case than in the FSRQs 
one. This is because the former sample is rather local z ^ 1 and 
probes an epoch in which the cosmological constant drives the ac- 
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Figure 13. CCFs estimated from the WMAP7 ILC map and the 21 -month 
Fermi-LAT EGB map with \b\ > 20° in three energy bands. The three sym- 
bols refer to 3 energy cuts E > I GeV, £ > 3 GeV (upper panel) and i? > 30 
GeV (below panel). Model predictions for different types of sources are rep- 
resented by continuous curves: FSRQs (black, continuous), BLLacs (red, 
dashed), star-forming galaxies (blue, dot-dashed). 
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Figure 14. Angular cross-power spectra of the EGB with the CMB. Dif- 
ferent hne-styles characteiize different models. Black, continuous curve: 
FSRQs. Red, dashed: BLLacs. Blue, dot-dashed: star-forming galaxies. 



celerated expansion in which the late ISW effect sets in. FSRQs, in- 
stead, preferentially sample a high redshift, matter-dominated, flat 
universe from which we expect only a very weak ISW signature. 
Once again, the difference between the BLLacs and star-forming 
galaxies models reflects the different biases of the two populations. 
The difference among the models is best seen through their angular 
cross-spectra, shown in Fig. [14] As for the auto-power case, we no- 
tice that the FSRQs cross-power peaks at smaller angles than that 
of the other models, reflecting the different redshift sampled by the 
different classes of objects. The smaller power of the FSRQs model 
reflects the intrinsic weakness of the ISW signal, as anticipated. 
The drop at / ~ 100 is due to the late-time (low £) nature of the 
ISW effect. 



5.3 Cross-correlation with QSOs 

In Fig.[T5]we plot the estimated CCF between the Fermi-LAT EGB 
and the SDSS-DR6 QSOs maps compared with model predictions. 
Symbols and line-styles are the same as in Fig. [T3] The expected 
CCF amplitude is rather weak. This is due to the quite sharp, high 
redshift peaks of the quasars' dN/dz which preferentially picks up 
the correlation signal from high redshifts. With a large z value, the 
(1 + z)^^H^\z) dimming term in Eq. |9] weakens the coiTelation 
signal. The fact that the correlation signal comes from large redshift 
irrespective of the EGB model explored, is more clearly seen in 
Fig.[T6]which shows that the angular power has indeed been driven 
toward small angular scales. We notice here that it is still physical 
to consider, for example, the star-forming galaxies model in rela- 
tion to the cross-correlation with QSOs. Indeed, a cross-coiTelation 
analysis alone cannot unambiguously identify EGB sources since 
all of them trace, with different biases, the underlying Dark Matter 
field. So, there may be a positive cross-correlation with the SDSS- 
DR6 QSOs catalog, which hardly consists of star forming galaxies, 
even if the underlying sources of the EGB are star-forming galax- 
ies, since both of them trace the the Dark Matter field. 

All CCFs are consistent with zero. The weak (~ la) positive 
correlation in the innermost bin for the £ > 3 GeV case (which 
would be consistent with theoretical expectations), disappears for 
the different choices of the cleaning method and/or when increasing 
the \b\ cut. 

Note that, although the SDSS DR6 QSO catalog has very high 
efficiency in the selection algorithm, stars are point-like sources 
that inevitably contaminate the catalog. To compute the stellar con- 
tamination, we extract a large number (~ 8 x 10^) of stars from the 
SDSS DR6 survey in the magnitude range 16.9 < g < 17.1 using 
the CasJobsQ website and compute the CCF between stars and the 
Fermi maps c''*(6). We have carefully checked that the contribution 
from contaminating stars to the CCF between QSOs and Fermi- 
LAT maps is fully consistent with zero and can safely be neglected. 

5.4 Cross-correlation with NVSS galaxies 

The dN/dz of NVSS galaxies peaks at lower redshifts than the QSO 
ones. For this reason the expected CCF signal is larger than for 
QSOs (see Fig.ll7t and the cross-power shifts at larger angles (see 
Fig.1181). FSRQs have a weaker cross-correlation signal (and cross- 
power) because there is little overlap between the NVSS dN/dz and 
the redshift distribution of the y-ray emission signal (see Fig.[T}. 

The estimated CCFs are consistent with zero and with theo- 
retical predictions, despite the fact that all Fermi sources, including 
normal galaxies, are also NVSS sources. At the moment, this re- 
sults is not worrisome given that the expected signal still fits within 
the large eiTor bars. It will be indeed interesting to see if it persists 
with larger statistics. 

5.5 Cross-correlation with LRGs galaxies 

Fig. [19] shows that the expected amplitude of the CCF with the 
LRGs for BLLacs and star-forming galaxies is larger than that 
of NVSS galaxies because the LRGs distribution peaks at similar 
scales but is much sharper. As a consequence the peak of the angu- 
lar cross-power for BLLacs and star-forming galaxies is quite sim- 
ilar to the NVSS one with some more power at small angles (see 

' http://casjobs.sdss.org/dr7/en/ 
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Figure 15. CCFs estimated from the SDSS DR6 QSO map and the 21- 
month Fermi-LAT EGB map with \b\ > 20° in three energy bands. The 
three symbols refer to 3 energy cuts E > I GeV, E > 3 GeV (upper panel) 
and £■ > 30 GeV (below panel). Model predictions for different types of 
sources are represented by continuous curves: FSRQs (black, continuous), 
BLLacs (red, dashed) star-forming galaxies (blue, dot-dashed). 



Figure 17. CCFs estimated from the NVSS galaxies map and the 21-month 
Fermi-LAT EGB map with |fo| > 20° in three energy bands. The three sym- 
bols refer to 3 energy cuts E > I GeV, E > 3 GeV (upper panel) and £" > 30 
GeV (below panel). Model predictions for different types of sources are rep- 
resented by continuous curves: FSRQs (black, continuous), BLLacs (red, 
dashed) star-forming galaxies (blue, dot-dashed). 
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Figure 16. Angular cross-power spectra of the EGB with the SDSS-DR6 
QSOs. Different line-styles characteiize different models. Black, continu- 
ous curve: FSRQs. Red, dashed: BLLacs. Blue, dot-dashed: star-forming 
galaxies. 
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Figure 18. Angular cross-power spectra of the EGB with the NVSS galax- 
ies. Different line-styles characterize different models. Black, continuous 
curve: FSRQs. Red, dashed: BLLacs. Blue, dot-dashed: star-forming galax- 
ies. 



Fig. |20t . There is no curve for the FSRQs case since this model 
predicts zero cross-correlation signal because there is no overlap 
between the dN/dz distribution of the LRGs and the predicted red- 
shift distribution of y-ray signal. 

The measured CCFs is consistent with zero for E > 3 and 
£ > 30 GeV. With the lowest energy cut, £ > 1 GeV, we detect 
a positive correlation signal at 9 < 2° at ~ 2a confidence. This sig- 
nal is remarkably robust to cleaning procedures and Galactic cuts. 
This feature is also robust to the choice of the y-ray events since it 
is also present when we only consider the so-called front y-ray 
photons which have a significantly better PSF. 

Theoretical predictions agree with this signal at the ~ 1.5o 
level. The fact that is only seen at low energies may indicates that 



the sources that contribute to the EGB at low energies are brighter 
than expected in our power-law model, i.e. the bias of the EGB 
sources {by in Eq|9]l is larger than expected. More intriguingly, 
this may indicate that there is a transition in energy in the sources 
contributing to the EGB from e.g. galaxies or BLLacs at low en- 
ergy (which cross-correlate with LRGs) to FSRQs at high energies 
(which do not). 

A further alternative is that this signal may come from the 
cross-correlation with the y-rays contributed by sources coincident 
with LRGs but that are still too faint to be detected. Luckily, more 
statistics and better understanding of the diffuse foregrounds will 
help in the near future to better characterize this feature (for exam- 
ple using finer angular bins or lower energy photons) and check the 
validity of the previous hypotheses. 
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Figure 19. CCFs estimated from the LRGs map and the 21 -month Fermi- 
LAT EGB map with \b\ > 20° in three energy bands. The three symbols 
refer to 3 energy cuts E > \ GeV, £ > 3 GeV (upper panel) and £ > 30 
GeV (below panel). Model predictions for different types of sources are rep- 
resented by continuous curves: FSRQs (black, continuous), BBLacs (red, 
dashed), star-forming galaxies (blue, dot-dashed). 



Figure 21. CCFs estimated from the 2MASS map and the 21 -month Fermi- 
LAT EGB map with \b\ > 20° in three energy bands. The three symbols 
refer to 3 energy cuts £ > 1 GeV, £ > 3 GeV (upper panel) and £ > 30 
GeV (below panel). Model predictions for different types of sources are rep- 
resented by continuous curves: FSRQs (black, continuous), BLLacs (red, 
dashed) star-forming galaxies (blue, dot-dashed). 
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Figure 20. Angular cross-power spectra of the EGB with the LRGs. Dif- 
ferent hne-styles characterize different models. Black, continuous curve: 
BLLacs. Blue, dot-dashed: star-forming galaxies. 
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Figure 22. Angular cross-power spectra of the EGB with the 2MASS galax- 
ies. Different line-styles characterize different models. Black, continuous 
curve: BBLacs. Blue, dot-dashed: star-forming galaxies. 



5.6 Cross-correlation with 2MASS galaxies 

The result of the cross-correlation between Fermi-LAT EGB maps 
and 2MASS catalog (that represents the most local of our samples) 
confirms the trend of the other cross-correlation analyses: increase 
of the expected CCF amplitude predicted by the models of BLLacs 
and star-forming galaxies (see Fig. I21t . angular power that shifts 
toward larger angles (see Fig. 1221 and zero correlation expected for 
an EGB solely contributed by FSRQs. We note that, as expected, 
the theoretical angula r cross-spectrum i s in go od agreement with 
the one computed by I Ando & Pavlidoul ( |2009|) . The small differ- 
ences likely arise from the fact that in our estimate we did not ac- 
count for the angular resolution of the instrument and did not filter 
the angular power spectrum accordingly. 



The measured cross-correlation signal is consistent with zero 
at all but small angular separations and for £ > 1 GeV, where we 
detect a hint of positive correlation. However, the reality of this 
correlation signal is questionable for two reasons. On one side, we 
found that this signal is rather sensitive to the cleaning procedure 
and to the Galactic mask adopted. Also in the LRGs case we made 
a further check using front events only, but the sensitivity to data 
cleaning technique still persists. On the other side, this signal could 
be related to some possible systematic errors in the treatment of the 
2MASS catalogue which h as been advocated to settle some contro- 
versy in the ISW detection jFrancis & Peacoc5l201(]|) . 
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Table 2. test of the observed data against the ACDM and SCDM models. 



Sources 


Current Data (x^) 


Forecast (AX") 




ACDM 


SCDM 


Five years 


Ten years 


FSRQs 


19.94 


20.33 


0.7 


1.3 


BLLacs 


19.21 


19.45 


0.6 


1.2 


Star-form. 


20.80 


21.50 


1.8 


3.5 



6 ANALYSIS 

In this section we quantify the capability of the coiTelation analy- 
ses presented in the previous sections to constrain the cosmological 
model and the nature of the sources that contribute to the EGB us- 
ing the observed ACFs and CCFs. For this purpose we only con- 
sider £ > 3 GeV maps that allowed to obtain robust results (in 
contrast with the E > \ GeV case) with reasonable errorbars (in 
contrast with the £ > 30 GeV case). 

To do this we have computed the following function: 

where rf^'* represents the measured CCF between the diffuse EGB 
and the catalog k (coinciding with the ACFs when k = j) and f^'*^ 
is the theoretical prediction. The i,j indexes run over the 5 angular 
bins of each of the various ACF/CCFs, while the k,l indexes run 
over the different ACF/CCFs. Ci,j,kj indicates the general covari- 
ance matrix obtained from the jackknife resampling, and contains 
the correlation between different angular bins as well as between 
the various ACF/CCFs relative to different catalogues. In practice, 
this statistics compares the measured ACE and CCFs presented 
in Section[5]for our 3 model predictions (BLLacs, FSRQs and star- 
forming galaxies), considering all angular separations 6 ^ 10° and 
taking into account the covariance among the different ACE/CCE 
estimates. 

The purpose of this comparison is twofold. First, we want 
to check whether this analysis is sensitive to the presence of a 
cosmological constant (mainly through the EGB - CMB compar- 
ison). Second, within a given cosmological framework, we want to 
check our ability of discriminating among competing models for 
the EGB. To answer these questions we have performed two differ- 
ent tests. 

The results of the first test are summarized in Table |2] 
Columns 2 and 3 list the %- values obtained when we consider 
the "concordance" ACDM model adopted in this paper (second 
column) and the y} values obtained when we consider a CDM 
Einstein-de Sitter model with fini = 1.0 (SCDM, third column). 
The different rows refer to the three EGB models considered (col- 
umn 1). The X'^ values obtained for the two cosmological models 
are rather small. This similarity reflects the insensitivity of our cor- 
relation analysis to the presence of a cosmological constant term. 
The situation will improve significantly with the future Fermi-LAT 
data. In columns 4 and 5 we list the increase of the y}, relative to 
the values in column 2, expected after f^bs = 5 ™d fobs = 10 years 
of observations, respectively. This forecast has been obtained by 
assuming Poisson errors i.e. by scaling the 21-month errorbars by 
\/7(,bs- One can see that with fobs = 10 years one expects to dis- 
criminate a ACDM model form a SCDM scenario at about 2<7 con- 
fidence level if the EGB were mainly contributed by star-forming 
galaxies. This estimate assumes Gaussian statistics and refers to the 



case of one free parameter only: the value of Q.^. In fact the situa- 
tion is likely to be more favorable since future data will also allow 
to improve the Galactic foreground model. In this way we will be 
able to extend the congelation analyses to lower energy bands. In 
addition, one can improve the effectiveness of the statistics by 
carefully selecting the range of angular scales to be considered or 
by restricting the analysis to a few CCFs, among which the CCF 
with the CMB will play a crucial role since, as we have pointed 
out, the ISW signal is very sensitive to the underlying cosmologi- 
cal model. 

As a second test we have computed the y} varying the fraction 
of objects that contribute to the EGB, /j, in the interval [0,1]. In this 
case we do not use the values listed in Table [T]but let the contribu- 
tion of each source type vary between 0% to 100%. As before, we 
do not explore the case of a mixed contribution from different types 
of sources. The sensitivity to /j comes, in this case, entirely from 
the linear dependence of the normalization of ACFs and CCFs on 
/j itself. We note that varying /j would affect the prediction for the 
y-ray flux as a function of redshift. We decided to neglect this effect 
given the large uncertainties in the modeling of the sources. The re- 
sults are displayed in Fig.|23] The plots show the y} as a function of 
the EGB fraction contributed by the source /j . The three panels re- 
fer to the three models explored: FSRQs, BLLacs and star-forming 
galaxies and, within each panel, the three curves refer to different 
observation times: 21 months (continuous, black), 5 years (dashed, 
red) and 10 years (dot-dashed, blue). The vertical lines are draw in 
con'espondence of the reference values of /j listed in Table [T] The 
minima of the y} do not coincide with the fiducial values for /j . 
However, with the current 21-month data the discrepancy is hardly 
significant (barely above l-o in the BLLacs model). We note that if 
one removes the physical constraint /j ^ 1 then for the FSRQs the 
minimum of the y} would be found at /j > 1 . The fact that FSRQs 
favor large values of /j is not surprising. It simply reflects the fact 
that their y-ray flux is preferentially produced at large redshifts (see 
Fig.[D and, as a consequence, that all their angular spectra plotted 
in Section |5] have smaller amplitudes that in the BLLacs and Star 
Forming galaxy cases. Notice, however, that overall the sensitiv- 
ity to the FSRQs component is very weak and it is not expected to 
provide significant constraints even after 10 years of data taking. 
The constraints will improve significantly with the future Fermi- 
LAT maps: in the case of BLLacs and star-forming galaxies, after 
10 years of data taking one would be able to reject the hypothesis 
of 100% contribution to the EGB at the > 3o confidence level (but 
only to ~ 1.5o for the FSRQs). 



7 DISCUSSION 

The results presented in Section|5]show that with the current uncer- 
tainties it is not possible to discriminate among the different models 
of source populations, considered in this paper, contributing to the 
EGB emission. The measured ACE and CCFs are generally consis- 
tent with zero and, in particular, the lack of a measurable ISW effect 
prevents us from infeiTing the presence of a cosmological constant. 
The y} test presented in Section |6l however, demonstrates that the 
situation will improve with the duration of the Fermi-LAT mission 
that, thanks to the sheer number of the collected photons, will allow 
to reduce the Poisson noise that contribute to the total uncertainties. 

However, it is worth pointing out that, at present, the main fac- 
tor limiting the correlation analyses is not photon counts statistics 
but the imperfect knowledge of the diffuse Galactic background 
which prevents us from using the data at very low Galactic lati- 
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Figure 23. Limits on the fraction of tliree possible EGB sources from tlie current 21-montli Fermi-LAT data (black, continuous) and the future 5-year (red, 
dashed), 10-year (blue, dot-dashed) Fermi-LAT observations. The vertical lines denote the maximum values of fractions from the cun'ent theoretical constraints. 



tude and, especially, in the lower energy band. Clues that on some 
scales the current Galactic diffuse emission is not perfectly mod- 
eled are found when we compute the ACF and when we cross- 
correlate the 2MASS galaxy and the QSO-SDSS catalogs with the 
residual Fermi-LAT maps obtained with the lowest energy thresh- 
old E > I GeV. In those cases we did find a weak correlation signal 
at small angular separation 6^2° that, however, is very sensitive 
to the cleaning procedure and that disappears when a progressively 
larger strip across the Galactic plane is excluded from the analysis. 
A more convincing coiTelation signal is detected at the same sep- 
arations and energy threshold when we cross-correlate the Fermi- 
LAT maps with the distribution of the SDSS LRGs. This signal is 
remarkably robust to masking and cleaning procedures. However 
this signal is above theoretical expectations and is only detected 
at low energy. If genuine, a positive cross correlation signal indi- 
cates that sources contributing to the EGB are traced by LRGs and 
prompts some further speculations. The fact that at £ '-^ IGeV the 
measured signal is larger than expected may indicates that, at this 
energy, the bias of the EGB sources {by in Eq|9) is larger than ex- 
pected. In addition, the fact that the cross-correlation signal drops 
for £ > 3 GeV may suggests that EGB is not contributed by a single 
type of objects. Instead, the EGB would be produced by different 
types of sources and their relative contributions depend on energy. 
In this particular case, the decrease at £ = 3 GeV could indicate that 
the relative contribution of FSRQs, whose cross-correlation with 
LRGs is expected to be zero, increase with energy. Another possi- 
bility is that the y-rays are physically associated to LRGs just below 
the Fermi-LAT detection threshold. 

The model of Galactic diffuse emission is expected to improve 
in the near term, when future dedicated analyses, based on a better 
photon statistics, will be available. As a result it will be possible to 
confirm or disprove these correlation signals and the validity of the 
above speculations. Moreover, improving the Galactic model will 
significantly increase photon statistics by extending the correlation 
analysis to lower energy bands. 

The actual reduction of the enor bars will be more severe 
than what expected from simple statistics since a larger sky cov- 
erage will allow to sample new structures, effectively reducing the 
Cosmic Variance. Future cross-correlation studies will al so bene- 
fit fro m plann ed galaxy reds hift surveys lik e EUCLID ([Laureiisl 
l2009h . JDEM ( lGehrelsli2O10l) and Big Boss jSchlegel et alj|2009l) 
that will be able to trace the large scale structure over a large frac- 
tion of the sky probing the crucial epoch in which the current ac- 
celerated expansion of the Universe has presumably started. 

The main aim of the tests presented in Section |6] was to 
illustrate the limitations of the correlation analysis of the current 
datasets and to show that future datasets will probably allow to 



marginally detect the EGB correlation signal and the ISW effect. 
Another way to improve the statistical significance of the future de- 
tections is to optimize the tools to compare model with data and to 
consider additional objects' catalogs in the cross-correlation analy- 
sis. Exploring the potential of the various statistical tools that could 
be employed in a cross-correlation analysis is beyond the scope 
of this paper. However, even with the paradigmatic case of the 
X~ statistics, significant improvement can be obtained by restrict- 
ing the comparison to some particular range of angular scales, en- 
ergy range or to some particular catalogs of objects. From a more 
general perspective, the most rewarding way of extracting astro- 
physical and cosmological information from the diffuse EGB is 
to perform a combined likelihood analysis combining the cross- 
correlation analyses presented in this work with the one-point 
statistics, angular auto-correlatio n and spectral information, as sug- 
gested bv lDodelson et al.l ( |2009|) . 

We stress that the EGB models used in this paper can and 
will improve significantly in the near future. Indeed, the number 
of resolved extragalactic sources is bound to increase with time. 
On one hand, this will allow us to resolve an increasing fraction of 
the diffuse EGB. On the other hand, the next generation catalogs 
of extragalactic source will allow to probe more reliably the faint 
end of the logN-logS of blazars and will provide better and bet- 
ter estimates of their contribution to the diffuse EGB. Population 
studies with higher statistics will also likely result in the detection 
of the expected break in the luminosity function, allowing to drop 
the assumption of a minimum Lmin- With strong constraints on the 
blazars' contribution, one could afford exploring the case of mixed 
contribution from different types of sources, especially if the anal- 
ysis can be performed over an extended range of energies, since 
the contribution from star-forming galaxy is expected to increase 
below I GeV. In this context one may also include the possible 
y-ray photons from DM annihilation in extragalactic halos, whose 
contribution to th e EGB has been modeled by seve ral authors (e.g. 
lUllio et all ( l2002h : lAndd ( l2005h : |Pieri et al.l l2009h ). 

Finally, it is worth pointing out that investigating the ISW 
effect using y-ray data (even if on its own it does not provide 
tight constraints on the cosmological parameters) is of consider- 
able importance since it would represent a consistency check for the 
ACDM model obtained from a new, independent dataset. Current 
ISW estimates that rely on different tracers for large scale structure 
are someti mes in tension w ith the amplitude value for a ACDM 
model (e.g. IHo et al The new analysis performed in this 

paper can contribute to quantify to what extent such measurements 
are reliable and identify possible systematic effects. 
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8 SUMMARY AND CONCLUSION 

In this work we have extracted the maps of the diffuse EGB from 
the 21 -month Fermi-LAT data by masking y-ray point sources from 
the Fermi-LAT data and by subtracting the various available models 
of the Galactic diffuse signal. 

These residual sky maps have been used to compute the an- 
gular two-point auto-correlation function of the diffuse signal. To 
minimize possible systematic effects that may affect the correla- 
tion analysis, we have thoroughly checked for the presence of spu- 
rious signals due to an imperfectly subtraction of the Galactic fore- 
ground. In practice, we have tested the robustness of all measured 
ACFs and CCFs to Galactic foreground models, cleaning proce- 
dures and masking strategies. In doing so we have optimized a 
strategy with combined cleaning and masking procedures that al- 
lows to reduce systematic uncertainties. We note that, as expected, 
the ACF is much more prone to systematic effects than the CCF 
which, by contrast, is remarkably stable. 

In addition to the auto-correlation analysis, complementary to 
the recent determination of the Fermi- LAT angular power spectrum 
bv lGomez Vargas et al.l ( I2OI0I) and bv lSiegal-Gaskins et al.l bOld) . 
we have cross-correlated these EGB maps with the WMAP7 CMB 
map, in an attempt to detect the ISW effect. Finally, to unveil the 
nature of the unresolved sources that may contribute to the EGB we 
have cross-correlated the Fermi-LAT EGB maps with the angular 
distribution of different types of extragalactic objects in several cat- 
alogs. More specifically we have considered the local population of 
2MASS galaxies, that should trace the EGB contribution from star- 
forming galaxies, and the more distant population of LRGs, NVSS 
radio galaxies and QSOs that may trace the population of blazars. 

All the measured ACFs and CCFs have been compared with 
theoretical predictions in which the mean EGB signal and its an- 
gular correlation properties have been modeled assuming only one 
type of contributing sources: star-forming galaxies, BLLacs or FS- 
RQs. The results of these comparisons allow, in principle, to con- 
strain the level of contribution to the EGB of the different sources 
and, mainly through the ISW signal, the value of fl^- 

The main results of our analysis can be summarized as fol- 
lows: 

• The measured ACF of the EGB above 3 GeV is consistent 
with zero at all angular separations. This result is in agreement, 
within the large error bars, with theoretical predictions since the 
auto-coiTelation signal is expected to be very small in all EGB mod- 
els explored in this work. In other words, the auto-cortelation func- 
tion does not seem to be the best statistical tool to reveal the nature 
of the sources that contribute to the EGB. In the E > \ GeV band 
we detect a ~ 2ct positive correlation signal at ^ 2° which, how- 
ever, disappears when we increase the size of the Galactic mask. 
Considering the sensitivity of the ACF analysis on the Galactic dif- 
fuse model, especially at low energy, we conclude that this signal 
is a spurious feature. 

• The measured ISW signal in the energy band is also consistent 
with zero at all energies. The expected amplitude, however, is sig- 
nificantly larger than that of the auto-correlation signal, suggesting 
that future data with better photon statistics and more accurate sub- 
traction of the Galactic contribution could allow one to detect the 
ISW signature. 

• The CCFs measured considering various objects catalogs are 
generally consistent with zero with a couple of exceptions: SDSS 
QSOs and 2MASS galaxies show a positive cross-correlation sig- 
nal with the £■ > 1 GeV Fermi-LAT photons for 9^2°. However, 
neither signal is robust to the cleaning procedure and the Galac- 



tic mask applied. A more intriguing correlation signal is found at 
these same separations and energy when we correlate the Fermi- 
LAT maps with the SDSS LRGs. Unlike the previous cases, this 
signal, which is detected at 2a, is remarkably robust to the model 
adopted for the Galactic diffuse signal, to the procedure to clean 
out spurious residuals and to the size and shape of the mask ap- 
plied to exclude the Galactic plane and the Bubble/Loop-I regions. 
The analysis of future Fermi-LAT data will clarify the reliability 
and the nature of all these features. 

• A simple test performed using all the measured correla- 
tion functions confirms, in a more quantitative way, that 21-month 
Fermi-LAT maps allow to neither investigate the nature of the EGB 
nor constrain the value of the cosmological constant. However, this 
analysis shows that 10-year Fermi-LAT data would allow to con- 
strain the contribution of star-forming galaxies of BLLacs to the 
EGB with a confidence level of > 3<J and to confirm the presence 
of a cosmological constant with a statistical significance of ~ 2 o. 

• These estimates are very conservative since they are based on 
a simple extrapolation of the 21 -month data assuming pure Poisson 
errors. The Galactic diffuse model, however, is likely expected to 
improve in the near future, while, new Fermi-LAT data will allow 
to better constrain the contribution to the EGB of some some ex- 
tragalactic objects like the FSRQs and the BLLacs. As a result, we 
will be able to extend the CCF analysis to lower energies, further 
improving the photon statistics, and providing more secure priors 
for the analysis. 
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